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Abstract 

The aim of this study is to compare the magnetic behavior of URu2Si2 under uniaxial stress 
along the a-axis with the behavior under hydrostatic pressure. Both are very similar, but uniaxial 
stress presents a critical stress a x smaller (0.33(5)GPa) than the hydrostatic critical pressure p x 
=0.5 GPa where the ground state switches from HO (hidden order) to AF (antiferromagnetic) 
ground state. From these critical values and from Larmor neutron diffraction (LND), we conclude 
that the magnetic properties are governed by the shortest U-U distance in the plane (a lattice 
parameter). Under stress, the orthorhombic unit cell stays centered. A key point shown by this 
study is the presence of a threshold for the uniaxial stress along the a-axis before the appearance of 
the large AF moment which indicates no-mixture of order parameter (OP) between the HO ground 
state and the AF one as under hydrostatic pressure. The two most intense longitudinal magnetic 
excitations at Qo=(l,0,0) and Qi=(0. 6,0,0) were measured in the HO state: the excitation at Qo 
decreases in energy while the excitation at Qi increases in energy with the uniaxial stress along 
the a-axis. The decrease of the energy of the excitation at Qo seems to indicate a critical energy 
gap value of 1.2(1) meV at <r". A similar value was derived from studies under hydrostatic pressure 
at p x . 

PACS numbers: 75.30.Mb, 75.50.Ec, 25.40.Dn, 78.70.Nx, 61.50.Ks, 75.25.Dk, 75.30.Kz 
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I. INTRODUCTION 



Puzzling heavy fermion physicists for more than twenty years, URu2Si2 is one of the 
most studied and least understood uranium compounds. The mysterious phase transition 
at To ~ 17.8 K of this 5/ heavy-electron compound is characterized by large bulk anomalies 
and sharp magnetic excitations in q-space and energy, at different Q- vectors (Qo=(l,0,0) and 
Qi = (0.6,0,0)). Concomitant with this order, a tiny but persistent antiferromagnetic moment 
(~ 0.02 hb) is measured in all samples with a wave-vector Q j4i r=(0,0,l) (equivalent to Q ). 
This tiny staggered moment is difficult to consider as the order parameter in a conventional 
antiferromagnetism frame as it cannot be reconciled with a jump of the specific heat AC /T ~ 
300 mJ/K 2 mol involving an entropy (S ~ 0.2Rln(2)). Because there is no determination of 
the order parameter (OP), the order in URu2Si2 is named the hidden order (HO). However, 
under pressure, URu2Si2 orders in a high moment antiferromagnetic (AF) structure with 
the wave-vector Qaf arid a moment of 0.36-0.4//^-. The well-defined phase diagram (see 
Fig{T]) shows that when URu2Si2 switches from HO to AF state at a critical pressure p x ~ 
0.5 GPa, the bulk superconductivity disappears^ as well the antiferromagnetic excitation 
Eo at Qo, signature of the HO phased At p^, the excitation Ei at Qi jumps from 5 meV to 
8 meV-. Under magnetic field (applied along c-axis), the pressurized AF phase is unstable 
and URu2Si2 reenters into the HO state-^(see Fig{T|). 

Thermal expansion measurements established that stress will increase T when it is ap- 
plied along the a-axis and decrease To when it is applied along the c-axis. The opposite 
effect is observed for the evolution of superconductivity*^. Recently a theoretical mode^ 
indicated that the space group of the PM (paramagnetic) state and of the HO state are 
different but may keep the atomic positions which explains why most of the local probes did 
not detect any modification in the crystal structure. This new point of view led us to study 
the variation of the crystallographic structure and revisit the magnetic behavior of URu2Si2 
under uniaxial stress along the a-axis. In a first experiment, Yokoyama et al. have shown 
that under uniaxial stress along the a-axis and below T the HO state is gradually mixed 
to the AF stated 

Uniaxial stress combined with hydrostatic pressure can bring information about the ex- 
change integrals between magnetic atoms. For example when an AF compound switches 
to PM state under pressure, a study under uniaxial stress can determine which parameter 
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FIG. 1. Schematic phase diagram (T,H,p) of URu2Si2- The data come from Ref. 
GPa corresponds to the critical pressure, p* ~ 1.4 to the pressure where URu2Si2 transits directly 
from PM state to AF state. The superconducting phase is not presented to simplify the phase 
diagram. 

(lattice parameter or ratio of lattice parameters) governs the magnetic behavior of the sys- 
tem. However uniaxial stress is a tool which is not used very often with neutron scattering 
as there is a high probability to break the sample and as large crystals with good ratio of 
the height by the diameter must be selected to realize a homogeneous uniaxial stress condi- 
tions. In spite of this, neutron users prefer to use thin samples as for the first experiment 
on URx^Si^ with the difficulty of poor homogeneity. 

This paper is organized as follows. In section |TT1 we present the experimental set-up. The 
results are presented in section IIIII Section [TV] is dedicated to the discussion of our results 
and their comparisons with previous data. Finally concluding remarks are given in section 
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II. EXPERIMENTAL SET-UP 

Uniaxial stress was applied along the a- axis of two single crystals of URu2Si2 coming 
from different batches. The first one is a small sample with a perfect cylindric shape of 
diameter (d) 3.78 mm and height (h) 1.66 mm and the second one has a parallelepiped 
shape of surface 12 mm 2 with a vertical a-axis of 8 mm length. An important parameter to 



perform reliable uniaxial stress experiments is, as we will see latter, the ratio of the height 
by the diameter of the sample that we define as k = h/d: k ~ 0.5 for the small sample 
and k ~ 2.0 for the large one (In this case an average diameter was calculated). The small 
sample comes from the same batch as samples previously used for high fielcU^ and pressure 
measurements^. The second sample is from a new crystal grown in a tetra-arc furnace and 
annealed for 5 days at 1075°C. The samples were installed successively at the bottom of 
the Institut-Laue-Langevin (ILL) uniaxial stress stick between two foils of cadmium and 
gold to flatten the surface defects of the sample and loading platforms, and to reduce the 
friction when the sample is pressed. The stick was installed in an ILL-orange cryostat on 
the cold-triple-axis IN12, CRG spectrometer at ILL. 

The small sample was principally used to determine the nuclear structure with neutrons 
at k/=1.48A _1 . Monochromator and analyzer were put flat vertically and horizontally, 
in order to enhance the effect of the graphite filter. The collimators were: open-60'-60'- 
60'. Two filters: one Be filter on kj and one graphite filter on kj (the graphite filter was 
oriented such that the neutrons of wavelength A/2 were diffracted out by the reflection 
(006) of the graphite; the neutrons of wavelength A are not diffracted out as the reflection 
(003) of the graphite does not exist). The large sample was used for inelastic and elastic 
magnetic scattering measurements with ky =1.5A _1 . Vertically curved monochromator and 
horizontally curved analyzer were used with collimators: open-60'-open-open, one Be filter 
on k/-, plus a not-cooled Be filter which was installed on kj when magnetic Bragg peaks were 
measured. 

A complementary Neutron Larmor diffraction (NLD) experiment was performed on IN22 
(CEA CRG-beam-line at ILL) on the large sample to obtain the distribution and the tem- 
perature dependence of the lattice parameters a and c. NLD exploits the Larmor precession 
of the neutron spin within well-defined magnetic field regions to measure the particle's 
wavelength with high accuracy (Larmor encoding). The beam is initially polarized by reflec- 
tion on a Heusler-111 monochromator. The field is in practice simulated by pairs of radio 
frequency spin flippers (RFSFs) separated by a magnetically screened volume (Bootstrap 
technique described in Ref.jl^]). These devices consist of a rectangular coil producing a 
vertical static field Bq and containing another coil generating a smaller field oscillating in 
the horizontal plane at a frequency u r f = 7 n -Bo {in is the neutron gyromagnetic ratio). 
The length of the precession region is defined as four times the distance (bootstrap set-up) 
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FIG. 2. Sketch of a Neutron Larmor Diffraction setup. The scattering plane is horizontal. Radio 
frequency spin flippers (in blue) are rotated so that their faces are parallel to the lattice planes 
which generate the process of diffraction. This peculiar arrangement cancels out the effect of beam 
divergence in the scattering plane. 



between the first and second RFSF in each arm of the spectrometer, respectively Li and Lj 
(see Figj2]). Larmor encoding is used to perform high resolution diffraction measurement 
thanks to the symmetry of Bragg's law which selects the neutron's wavelength for a par- 
ticular lattice spacing d. Consequently, the total phase, namely the rotation angle of the 
magnetic moment of the neutron, is denned as tp = u r f 4 (Lj + Lf) m n /h ■ d ■ sin(6 , s) and 
measured by projection on the quantization axis of a Heusler analyzer (m n is the neutron 
mass and h is the Planck's constant). Rotating the RFSFs so that their faces are parallel 
to the lattice planes generating diffraction ensures that the latter relation is fulfilled for any 
wavelength in the incoming beam's bandwidth. One finally obtains the very simple equality 
Sip/ifo = Sd/do where subscripted variables are reference quantities. The advantage of this 
method is that the high resolution is achieved with modest beam collimation as a d-spacing 
creates the same phase for all neutrons. The experiment was performed with an incident 
wavelength kj=2.662 A -1 . Using a Larmor frequency u r f = 670 kHz and a total effective 
length (4(Lj + Lf)) of 3.44 m, the total Larmor phases p was close to 8650 rad with a 
theoretical resolution of the order of 5 ■ 10~ 6 . The technical configuration is described in 



Ref. 
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III. RESULTS 



A. Elastic nuclear scattering 

URii2Si2 crystallizes in the body-centered-tetragonal structure (space group: I A/ mmm). 
Only 3 parameters are needed to describe the nuclear cell: a and c, the lattice constants and 
zsi, the silicon atomic position along the c-axis. It is now believed^ that the symmetry goes 
from body-centered-tetragonal in the paramagnetic state to simple-tetragonal in HO. Both 
space groups have the particularity to keep the same atomic positions so to conserve the 
same nuclear structure factor in both states. As uniaxial stress is applied along one a-axis 
(named a^, v for vertical), the tetragonal symmetry is broken and the second a-axis (named 
a.h, h for horizontal) is no longer equivalent to a^: the symmetry becomes orthorhombic. As 
measurements are performed on a triple-axis spectrometer only few nuclear reflections are 
available in the plane (si h ,c). However, the nuclear reflections Q =(1,0,0) and Q=(0,0,1) 
stay forbidden under uniaxial stress in the paramagnetic state (above T ) which indicates 
that the translation (1/2, 1/2, 1/2) is conserved in the crystal. To confirm this assumption, 
it must be checked that the nuclear reflection Qo=(0,l,0) stays forbidden. This means that 
only the 4-fold axis disappears and the space group becomes / mmm in the paramagnetic 
state under uniaxial stress. In this new space group the nuclear structure is described with 
a new parameter: the position along c-axis of the ruthenium atom (z# u ). 

It is interesting to note that, as the space group 1 4 /mmm, the space group of the PM 
state (Immm) under uniaxial stress along the a-axis (a a ) may also be transformed in the 
HO state into a space group with the same atomic positions. The space group of the HO 
state would then be Pnnm or Pnnn. This may indicate an important property of the U 
site symmetry: the transition into HO state does not change the atomic position and HO 
state exists even if the four-fold axis is broken to become a two-fold axis. 

The pressure dependence of the lattice parameters were determined at 8K on the small 
sample using the positions of the nuclear Bragg reflections Q=(0,0,2) and Q=(1,0,1). Both 
lattice parameters in the scattering plane increase slightly and almost linearly with the 
pressure as expected in the elastic deformation regime. As the increase of dc/da a is ap- 
proximatively 2.5 larger than dah/da a which is almost the ratio c/a, the ratio c/a^ remains 
almost constant under uniaxial stress for the range of pressure from 0.2 to 0.65 GPa. Com- 
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pared with the dependences of the cell parameters deduced from the elastic constants (given 



in Ref. 



16|), there is a quite good agreement. 



The intensities of the nuclear reflections at Q=(0,0,2) and at Q=(1,0,1) increase by 77% 
and 88% respectively at the maximum of the applied uniaxial stress. Even if these reflections 
have a small structure factor compared to the nuclear structure factor of the largest intense 
reflection Q=(2,0,0) (1/5 and 1/12 respectively), they are largely affected by the extinction 
because of the size of the sample and the long wavelength of the neutrons. However the 
larger increase of the weaker reflection indicates a modification of the nuclear structure 
factor certainly due to a modification of zsi- This effect will be studied in a future dedicated 
diffraction experiment. 



B. Elastic magnetic scattering 

The magnetic Bragg peak measured at Qo and at 30 K, on the large sample, presents no 
intensity for the whole range of measured stress. Transverse scans are presented in Figj3] for 
different stress at low temperature and at T=30K for the largest stress we applied: <r a =0.3 
GPa. The resolution of the spectrometer of 9.12(2)-10 _3 A~ 1 was determined using the widths 
of the two nuclear peaks at Q=(1,0,1) and Q=(0,0,2). To determine the magnetic correlation 
lengths, the magnetic peaks were fitted as a convolution of the gaussian resolution with a 
lorentzian function. The magnetic intensities given by neutron scattering at ambient pressure 
can be associated to the AF volume fraction assuming that it comes from residual AF 
components induced by local defects inducing local stresses. At low temperature (T=2K), 
the magnetic intensity increases slightly for stresses lower than 0.2 GPa and drastically for 
larger stresses. The stress dependence of the AF volume, assuming a saturated AF moment 
in the AF state of 0.36±0.04 /xg (mean value of AF moment of URu2Si2 in the AF phase 
under hydrostatic pressure), and the magnetic correlation length at low temperature are 
presented in FigJH Both start to diverge at a a = 0.3 GPa and we estimate the critical stress 
a°j. ~ 0.33GPa which is defined as half of the crystal volume in the AF state. With the 
hypothesis of a linear relation between magnetic intensity and AF volume, the repartition 
of volumes is 30%/70% of AF/HO states at cr a =0.3 GPa. A critical temperature T£ can 
be defined as well when half of the crystal volume in the AF state: T£ will arise from 
K just above a%- However, this estimation is based on the assumption of identical AF 
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moment values when uniaxial or hydrostatic pressure is applied, which may be wrong as the 
tetragonal symmetry is broken with the stress. 

The temperature dependence of the AF moment presents different behaviors according 
to whether the state is HO or AF. This modification happens at a value of the AF moment 
around 0.06 \xb and can be chosen as a criterion for the transition temperature T^ 7 from 
HO to AF order (see inset figure [3]). T£ is larger than as it corresponds to a magnetic 
intensity equivalent to 0.25 fis- Similar behavior was found in previous measurements under 
hydrostatic pressure^^^. 

A larger critical stress ~ 0.55 GPa was determined for the small sample. This difference 
is due to the reduced free expansion condition when the ratio k is small. It is well-known 
that for samples with small k ratio (n < 2) as for our small crystal or for samples of M. 
Yokoyama et alM (n ~ 0.2), the experimental conditions are between uniaxial stress (free 
lateral expansion) and uniaxial strain (no free lateral expansion) which, in the last condition, 
increases the critical pressure. The behavior of the AF volume versus uniaxial stress in the 
sample of the Ref. 12j (plotted to comparison in Fig Hj) cannot only be explained by the 
experimental conditions. The increase of the AF volume is largely affected by the quality 
of the sample: it shows a large AF volume already at low hydrostatic pressure^ which does 
not correspond to the usual results^. 



C. Inelastic magnetic scattering at Qo and Qj nc 

For the inelastic neutron scattering measurements at Q , the variation with stress of 
the intensity of the excitation is small (see Figj5]) and the energy gap Eo decreases lin- 
early from 1.68(l)meV down to 1.27(l)meV. To evaluate the variation of this intensity some 
assumptions were made. First, we consider that the intensity of the excitation is related 
to the volume of the HO phase, which means that this excitation is characteristic of the 
HO. Secondly, URu2Si2 can be described by a single- mode approximation: sharp dispersion 
at the antiferromagnetic position Q , which gives the intensity of the magnetic excitation 
~ I/E'er^. Then the variation of the intensity as a function of uniaxial stress should be 
I(cr a )=I(cr a =0)V HoEo(a a =0) / Eo(a a ) . This gives only a slight decrease of the intensity in 
agreement with the inelastic neutron scattering measurements. There is no large modifica- 
tion of the intensity because at the larger stress we applied most of the sample was still in 
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FIG. 3. Intensity in log scale of the magnetic Bragg peak at Qo for different uniaxial stress 
at low temperature T=4K and above To(T=30K) for the largest stress. The inset displays the 
temperature variation of the magnetic intensity (also in a log scale) and shows the different regime 
according URu2Si2 is in the HO or in the AF states. T£ is defined in the text. 
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FIG. 4. Uniaxial stress dependence of the AF volume and of the correlation length for our two 
crystals and also the AF volume of the crystal determined by M. Yokoyama in Ref . . The AF 
volume is determined with the assumption that the moment in the AF state is m.Ap>=0.36fiB- For 
indication, the correspondence between AF volume and an intrinsic moment value is given, but 
noticed that this scale (in red) is not linear. LS and sS mean large and small samples respectively. 



the HO state. The evolution of the energy gap and of the intensity versus uniaxial stress 
along a-axis are very similar to the results under hydrostatic pressure^. Considering the 
critical stress <r^=0.33 GPa, determined from the elastic magnetic measurements, we can 
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FIG. 5. Magnetic excitations at Q =(1,0,0) and Qi=(0.6,0,0) for a a from to 0.3 GPa at low- 
temperature. 

deduce a critical energy gap Eo_ C rit of 1.2 meV at a% which has the same value than the 
critical energy gap under hydrostatic pressure at p x . However the critical pressure p x is two 
times larger (0.5-0.7 GPa£>^) than the critical stress a®. 

For the incommensurate excitation at Qi, the energy gap slightly increases with uniaxial 
stress then its intensity slightly decreases (see Figj5]). The increase of the energy gap from 
4.27(1) meV to 4.48(1) meV is again similar to the results with hydrostatic pressure^. Using 
the single-mode approximation, a larger decrease of intensity is expected. However, contrary 
to the excitation at Qo, the incommensurate excitation at Qi does not vanish in the AF 
state: its energy gap Ei just shifts to higher energy. Therefore as the HO-AF transition is 
first order, it is possible that we are measuring a mixture of E x in the HO state at an energy 
transfer of ~4.5meV, with Ei in the AF state at higher energy. This may explain why the 
intensity decreases less than expected. In the case of hydrostatic pressure, the energy gap 
Ei jumps from ~ 5meV to ~ 8meV— , but it may be smaller in the uniaxial stress case. 
Nevertheless, the initial evolutions of the gaps under uniaxial stress or hydrostatic pressure 
below the stress or pressure threshold in the HO state, are very similar and give the same 
critical energy gap E _ crit ~1.2 meV. 
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D. Thermal expansion and distribution of the lattice parameters a and c 



Using NLD, we have measured the temperature variation of the a and c lattice parameters 
in URu2Si2 at ambient pressure. Thermal expansions from 10K to ~ 85K along the two 
axis have been inferred from the (2,0,0) and (0,0,4) Bragg reflections respectively. They 
were performed with increasing temperature. The thermal expansion value at 10 K given 
by dilatometry-22 was used to normalize our data. The results are shown on Figj6] and 
compared to a dataset obtained by means of a three-terminal capacitance method (from 



22 and 



231]) . Along the a-axis, we observe two clear changes in the slope at T ~ 17.5 K 
and T x ~ 45 K. The first one corresponds to the transition from HO to the large moment 
AF phase. The second one is located at the position of the maximum in the easy-axis 
magnetic susceptibility^ 4 -. Along c, To is marked by an inflection point in the curve whereas 
T x corresponds to a local minimum. The two diffraction and bulk measurement datasets 
are only qualitatively consistent: the main features are seen but absolute values differ. We 
have no clear explanation for this discrepancy, which might come from a different sample 
thermalization. 

Another consequence of the linear relation between the Larmor phase ip and the lattice 
parameter d is that the final beam polarization, Vj((fo), yields the cosine Fourier transform 
of the lattice spacing distribution function f(d/d ) through : 

(-«)> 

where f(d/do) is assumed to be gaussian, V% is the beam polarization extrapolated to zero 
phase, ipo is the total phase corresponding to the mean neutron wavelength (ip Q oc <i oc A) 
and Sd/do is the FWHM of f(d/do). In practice, we find that Vi ~ 0.72 while the natural 
polarization provided by reflection on the Heusler monochromator and analyzer is ~ 0.90. 
The difference between usually natural polarization and Vi is due to the divergence of the 
neutron beam through RF coils. 

Fig H shows the remaining normalized polarization (Vf(ipo)/Vi) for <£>o=8650 rad at T 
= 2 K. This polarization gives a relative distribution width 8a/ a$ = 4.5(3) • 10~ 4 and 
Sc/cq = 4.4(2) • 10~ 4 . A recent published experiment using NLD^, has determined that 
the distribution of the lattice parameter c was surprisingly two times smaller than the 
distribution of the in-plane parameters a, i.e. 2.1 ■ 10 -4 and 4.05 ■ 10~ 4 respectively. The 
outcome of our study is that both distributions are isotropic, independent of the temperature 
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FIG. 6. Thermal expansions of a and c lattice parameters in URu2Si2- Dots, full lines and dashed 
ines represent our NLD data and the results from dilatometry measurements taken from [2| and 
23j], respectively: red and blue colors are for the a and c lattice parameters. 



between 2 and ~ 80 K, within error bars (see inset Fig. [7]). We would like to point out that 
those values have to be taken as high limits because extrinsic effects may slightly depolarized 
the beam and that the quality of our sample of URu2Si2 is as good as good as perfect high 
quality silicon single crystal (shown comparison on FigJTJ). Nevertheless, the values of lattice 
parameter distributions are the same along a and c. 

IV. DISCUSSION 



Within the linear elastic deformation regime, the hydrostatic pressure and uniaxial 
stresses are coupled with the strains by the well-known stress tensor which can be repre- 
sented by the matrix—. The pressure variation of the strains of all the U-U distances 
in the crystallographic cell as well as the ratio t]=c/a. were calculated using the elastic 
constants obtained with the ultrasonic-sound velocity extrapolated to T = 0K-. Table U 
summarizes the derivative coefficients which are relevant for our study. 

Before going further, Table U can be only used in the case of hydrostatic pressure (column 
2) or uniaxial stress along a^ in the free lateral expansion conditions (column 3), as empha- 
sized in Ref. 2^] (study made in the case of isotropic crystal). The free lateral expansion 
condition is fulfilled only if the ratio k ^ 2. For our small crystal or any crystal with k < 2, 
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FIG. 7. Final beam normalized polarization as a function of the total phase (fQ. Black line 
represents the sake of comparison, a resolution curve which has been obtained by scattering on 
a "perfect" Si single crystal. Red and blue points correspond to the value of the normalized 
polarization measured at (po = 8650 radians and T~2K. Red and blue lines correspond to a law 
of the form of Eq. Q] applied to the respective cases of a and c lattice parameter. Inset: Final 
beam polarization at tpo = 8650 rad as a function of sample temperature. 
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TABLE I. Derivatives of the strain in the linear elastic regime of a„, ijh and r\ v versus X with 
X={p,o" a }. The strain is defined as the relative variation of a parameter. It is represented with a 
hat on the parameter (£ = AC/C = ^ttw^ ) • 
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Table d cannot be used. In the case of very thin samples with no free lateral expansion, the 
set-up fulfills the experimental conditions for a uniaxial strain where can be estimated 
to ~ — I/en = —3.9 GPa -1 . In this condition, the critical stress for URu2Si2will increase 
by 12%. However as the experimental conditions are not under control if k <C 2, the results 
are largely uncertain. In the experiments of Ref. 121 ] . as n ~ 0.2 their results cannot bring 
any reliable output for the critical stress. 

Let us now compare the hydrostatic critical pressure p x generally found in the range 
of 0.5-0.7 GPa to the uniaxial critical stress along a-axis a" ~0.33 GPa. At the HO-AF 
transition, we consider that the relevant parameter (() is either the U-U distance (a) or 
the ratio t]. At the hydrostatic critical pressure or at the uniaxial critical stress, ( should 
have the same critical value (( x )- This is also true for £ which is the relative variation of ( 
(C = A(/(= - ^ > ( ^°' > ) • So ( x = §^X X is a constant independent of constraint X (see figure 
in Table HIT) where X may be either hydrostatic pressure (p) or uniaxial stress (a a ) (X={p or 
<r a }). To fullfill this condition, at the transition when URu2Si2 switches from the HO to the 
AF state, the experimental ratio ~ 1.7±0.3 should match the ratio J^-/§^ calculated in 
Table [Til Only a,; fulfills this relation and the shortest U-U distance in the plane appears as 
the best candidate to control the magnetic properties of URu2Si2- However to complete our 
discussion, the parameters rjh and r) v have been also considered in the following. In Table [TTJ 
( x have been calculated, at the critical hydrostatic pressure p x and at the critical uniaxial 
stress a%. It is interesting to note that with the shortest U-U distance in the plane as the 
relevant parameter, if the uniaxial stress is applied along the [1,1,0] direction, the critical 
stress a% x should be larger 0.8 GPa) than the hydrostatic critical pressure or uniaxial 
critical stress along a-axis. Also it is not possible to switch from HO state to AF state 
applying a stress along the c-axis in agreement with thermal expansion results. 

The NLD result shows a difference between our distribution of the lattice parameter 
c and the previous distribution measurement then there is almost no difference for the 
distributions of the lattice parameter a^. Their c distribution is two times smaller than 
their a distribution, whereas we obtain an isotropic distribution along these two directions. 
Nevertheless the tiny AF moment does not reveal a large difference, namely 0.020(4) /ig in 
our case to be compared to 0.012 hb in their case. Our larger distribution may be explained 
as we use a larger crystal ~ 100mm 3 to realize reliable inelastic experiments. 

According the phase diagram under pressure, an intrinsic AF moment exists only in the 
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c 


Co 


fir 1 Fit 

da a 1 dp 


Cx at Px 


Cx at a% 


Vh 


2.318 


0.5 


0.72-1CT 3 


0.20-10" 3 


Vv 


2.318 


4.5 


0.72-10" 3 


1.78-10- 3 


a v 


4.125A 


1.6 


-1.68-10 -3 


-1.45-10 -3 




TABLE II. Table of the relevant parameters C = Vh, Vv and a v : value at ambient pressure, ratio of 
the derivatives of C between uniaxial stress along a-axis and hydrostatic pressure, and values of Cx 
at the critical pressure p x =0.6 GPa and at <r"=0.33 GPa. The figure on the right explains the 
assumption why Cx should be equal at p x and cr", the critical value where URu2Si2 switches from 
HO to AF. 



AF state (with a value hi^f — 0.36-0.40 fiB at low temperature^). This assumption means, 
in the regime of linear elastic deformation, that the AF moment exists only above (or below) 
a critical value Cx of the relevant parameter. Thus it is possible to calculate the AF volume 
which corresponds to the integration of the distribution of the relevant parameter above (or 
below) to the critical value Cx- 111 this same way, the variation under hydrostatic pressure 
or uniaxial stress of the AF volume can be calculated: the variation of C, which corresponds 
to the average value of the distribution is given by the formula: £(X) = Co(l + ^X), with 
X={p or a a }. With the assumption of a gaussian distribution of the relevant parameters C, 
the AF volume is given by: 



vwv =/ ^ ( - ] dC 

2y / hT2^(X g - X) 
5C/Co 




(2) 



where Vo is the total volume of the sample, and Co is the value of the relevant parameter C 
at ambient pressure. The erfc function is the complementary error function which can be 
found in the scientific library of Python^ (Numpy and Scipy). 

Formula [2] shows that the slope of AF volume versus the pressure or the uniaxial stress 
increases when width of the distribution dC is smaller but also is larger for uniaxial stress 
than for hydrostatic pressure as the derivative of the strain ^ is larger. This may explain 
why with samples of bad quality (large distribution 5Q the moment has a smaller slope and 
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increases linearly with pressure. 

FigJS] represents the pressure and uniaxial stress variation of the AF moment (m^ * 
a/VafG°) /Vq) with or 77 as relevant parameters. Our data under uniaxial stress along 



a-axis are compared to the hydrostatic pressure results of [25|. It is clear that 77 does not 
appear to be the relevant parameter. With a^ as relevant parameter, the shape of the 
transition at 0% or p^ is well explained. However, it is not possible to explain anymore the 
tiny moment at ambient pressure with only a simple gaussian distribution. We manage to 
reproduce the pressure dependence of the tiny AF moment with two other sources: either 
with an additional contribution to the a distribution (a tail in the distribution due to defects) 
or to intrinsic exotic nanostructured defects related to the high sensitivity of URu2Si2 to 
pressure, uniaxial stress or magnetic field^. 

The suspicion that these exotic defects come from the unusual effects is due to the fact 
that the detected tiny AF moment appears related to energy gap at Qo, which is a signature 

Eq(X) 

of the hidden order-, and to the AF moment hiaf by the relation: hio(X)=hiaf e ~ 
with E r ~ 0.64 meV and X corresponding either to magnetic field studies at p=0 or uniaxial 
stress at H=0 (Figj9]). The hydrostatic pressure dependence of the AF energy gap E was 
previously determined at low temperature^: the energy gap Eo(p) decreases linearly with 
pressure up to the critical pressure with a critical gap E _ cr i t ~ 1-1(1) meV. In this study, 
the uniaxial stress dependence of the AF energy gap E has the same behavior with a larger 
slope but a smaller critical stress 0% leading to the same critical gap E _ cr i t — 1-2(1) meV. 

A first order transition at cr£ or p^ indicates usually a strong repulsion between the two 
order parameters governing each side of this transition line where both order parameters can 
be mixed only if a coupling exists between them and in particular if both break time reversal 
symmetry. The invariance of the nuclear crystallographic structure between the PM and HO 
states with the loss of symmetric elements in the HO state keeping the atomic position at 
ambient pressure as well as with hydrostatic pressure and also under uniaxial stress along a, 



is in agreement with the idea developed in Ref. 
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and 



331 ] with the proposals for quadrupole 



or hexadecapole solutions for the hidden order. The ordering of any even-parity multipole 
will not break time reversal symmetry and thus at low stress mixing between HO and AF will 
not occur. However two resonant x-ray scattering measurements ruled out the possibility of 
any quadrupole ordering by resonant x-ray scattering^^. A hexadecapole remains a sound 
solution. Such a ground state was proposed in a model based on a unified complex order 
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FIG. 8. Pressure and uniaxial stress dependences of mo with rj or a. v as relevant parameters. The 
fits were made with imiF=0.4//B and using the formulaEl The dashed and full lines correspond to 
rj and a v as relevant parameters respectively. For the hydrostatic pressure fit, the parameters were 
25l ]. for uniaxial stress they come from this study. The derivatives were taken 



taken from Ref. 
from Tabled! 



parameter where the real part is a hexadecapole (not breaking time reversal symmetry) 
and the imaginary the magnetic dipole (breaking time reversal symmetry)^^. However, it 
was proposed that the critical stress along a (or [1,1,0]) should be two times smaller than 
the hydrostatic critical pressure^. It would be interesting to probe the critical stress along 
[1,1,0] direction, as we expect here a large increase of the critical stress compatible with the 
shortest U-U distance as the relevant parameter (~ 0.8 GPa). 

On the other hand, mixture of HO state with AF state is possible in models with odd 
multipole ordering. The octupole model developed by Kiss and Fazekas^ seems in dis- 
agreement with the results of this paper as they have considered a coexistence of both order 
parameters under uniaxial stress. The dotriacontapole also breaks time reversal symmetry-^. 
Both order parameter, HO and m z belong to the same group representation and then may 



generally be mixed. However according to the Ref.[39|, the tiny AF moment is not a real 
AF moment as proved by the latest NMR results^ where neither the large (0.3/^b) nor the 
small (0.02-0.04/ie) moments were detected, excluding a homogenous small static classical 
AF moment in URu2Si2- Then the small signal measured by neutron scattering is the non- 
collinear dotriacontapole itself which has an extremely short-range stray field. However, such 
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FIG. 9. Tiny AF moment at low temperature versus AF gap Eo in a log scale from measurements 
under magnetic fiel d 13131132 and this study under uniaxial stress. The green area is not taken into 
account in the fit because it corresponds to a mixture of both phases (HO and AF). The red curve 
corresponds to a fit with E r = 0.64 meV and m^F = 0.36 \jlb- 



a multipole may only have a finite cross-section at high momentum transfer— which seems 
not compatible with the small momentum transfer of the tiny moment. Moreover, there is 



no experimental trace of the dotriacontapole in the AF state as proposed in Ref . [39J . Other 
contradiction: if the U-U distance is the relevant parameter, then the critical value for a 
occurs for a variation of 1.5%, one order of magnitude larger than the present experimental 
results (see Table HT|) . 

In previous considerations, the 5/ localized character of the U atoms plays a key role. 
Other possibilities exist for the hidden order parameter corresponding to m(Qo) based on 
more itinerant models such as a dynamic order parameter (symmetry breaking by dynamical 
antiferromagnetic fluctuations of the hidden order )^—. In this case, the HO breaks time- 
reversal symmetry only for a short time. The theoretical Fermi Surface computed for the 
symmetry-broken state suggests that the body-centered translation vector is broken in the 
HO phase, in good agreement with the observation and the fact that the Fermi Surface of HO 
and AF are quite similar, i.e. Qho=Qaf, as seen in the experimental results^. Calculation 
of Fermi Surface in the orthorhombic structure symmetry, using the cell parameters induced 
by uniaxial stress, may ascertain this itinerant model. 

18 



V. CONCLUSION 



Comparing the uniaxial stress along a-axis to hydrostatic pressure measurements, both 
phase diagrams are quite equivalent with a critical pressure almost two times smaller in the 
case of the stress (0.33 GPa compare to 0.6 GPa). The magnetic properties of URu2Si2 
appear to be governed by the shortest U-U distance (a parameter) and not by the ratio 
7] = c/a. NLD results invalidate the simple model of large lattice parameter distribution to 
explain the tiny AF moment as extrinsic. The study of the hidden order state (more exactly 
of its fingerprint, the excitation E ) under uniaxial stress indicates that this order can exist 
in a four-fold axis local symmetry as well in a two-fold axis symmetry. It is in agreement 
with a loss of the four-fold axis symmetry on the U site when entering in the HO state as 
for the group PA 2 /mnm (This space group is one of two possible space groups proposed by 
H. Harimaii, but the only one which loses the four- fold axis on the U site). 

The most promising model is the hexadecapole model developed in Ref. 
further discussed using symmetry argument in Ref. 



36 



and [37J and 



331 ] . It will be interesting with this 



model to estimate the evolutions of the excitations at Qo and Qi and to compare them to 
experimental results already measured under pressure, stress and magnetic fielcr^i^. 
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